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The HELLAS2XMM Survey. XII. The infrared/sub-millimeter view 
of an X-ray selected Type 2 quasar at z « 2 
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ABSTRACT 

We present multi-wavelength observations (from optical to sub-millimeter, including Spitzer 
and SCUBA) of H2XMMJ 003357.2-120038 (also GD 158_19), an X-ray selected, luminous 
narrow-line (Type 2) quasar at z=1.957 selected from the HELLAS2XMM survey. Its broad- 
band properties can be reasonably well modeled assuming three components: a stellar compo- 
nent to account for the optical and near-IR emission, an AGN component (i.e., dust heated by 
an accreting active nucleus), dominant in the mid-IR, with an optical depth at 9.7 \im along 
the line of sight (close to the equatorial plane of the obscuring matter) of x(9.7) = 1 and a full 
covering angle of the reprocessing matter (torus) of 140 degrees, and a far-IR starburst com- 
ponent (i.e., dust heated by star formation) to reproduce the wide bump observed longward of 
70 i_im. The derived star-formation rate is w 1500 M Q yr . The overall modeling indicates 
that GD 158_19 is a high-redshift X-ray luminous, obscured quasar with coeval powerful AGN 
activity and intense star formation. It is probably caught before the process of expelling the 
obscuring gas has started, thus quenching the star formation. 
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1 INTRODUCTION 

Recent evidence that almost every local spheroid contains a super- 
massive black hole (SMBH) in its center whose mass is pro- 
portional to that of its host (e.g., Magorrian et al. 1998; Geb- 
hardt et al. 2000; Ferrarese & Merritt 2000) suggests that the 
growth of black holes and their host galaxies is mostly coeval. 
These observational results have found further support from hydro- 
dynamical simulations of galaxy formation, where active galac- 
tic nuclei (AGN) provide, mainly through winds and outflows, 
feedback between SMBHs and host spheroids (e.g., Cavaliere & 
Vittorini 2000; Di Matteo, Springel & Hernquist 2005; Menci et 
al. 2006, 2008 and references therein). In this context, the most 
likely scenario is one in which most of the massive galaxies pass 
through a luminous starburst phase (with star-formation rates SFR 
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of a few hundred Mq yr -1 ) before the SMBH has built sufficient 
mass (10 8 ~ 9 M for a lo 12 - 4-13 - 4 M Q dark-matter halo; Granato 
et al. 2004) to shine as a luminous quasar (e.g., Silk & Rees 1998; 
Fabian 1999; Hopkins et al. 2006). Initially the quasar would be ob- 
scured by the star formation that feeds it, but may eventually clear 
away gas and dust through various feedback mechanisms. Under 
this hypothesis, the fraction of obscured AGN can be interpreted as 
a measure of the timescale over which the nuclear feedback is at 
work (for a theoretical and observational investigation of this topic, 
see Menci et al. 2008; Fiore et al. 2009, and references therein). 
The SMBH emission will subsequently decline to leave a quiescent 
spheroidal galaxy. Sources caught during the initial luminous star- 
burst phase may comprise the bulk of the sub-millimeter (sub-mm) 
galaxies discovered by SCUBA over the last decade (e.g., Smail et 
al. 1997; Hughes et al. 1998), many of which being ultra-luminous 
infrared galaxies (ULIRGs; L > 10 12 L ) at z « 1-3 (e.g., Blain 
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et al. 2002); see also Sanders et al. (1988) for a similar scenario 
applied to ULIRGs and luminous quasars. 

Although the first investigations of AGN activity in sub-mm 
sources showed little overlap between the X-ray and the sub-mm 
source populations (e.g., Fabian et al. 2000; Severgnini et al. 2000; 
Almaini et al. 2003; Waskett et al. 2003), subsequent more sensi- 
tive X-ray observations in the Chandra Deep Field-North (CDF-N) 
were able to detect weak X-ray counterparts of sub-mm sources 
(e.g., Alexander et al. 2003, 2005a, 2005b). These studies, coupled 
with the availability of spectroscopic information (e.g., Chapman et 
al. 2003, 2005), revealed that the AGN contribution to the sub-mm 
emission is of the order of a few per cent, and that AGN in typ- 
ical sub-mm galaxies are growing almost continuously, although 
the SMBHs in these galaxies are possibly several times less mas- 
sive than seen in local massive galaxies (Alexander et al. 2008). 

Over the last few years, sub-mm evidence for the coeval 
growth of SMBHs and galaxy bulges has been reported by Page 
et al. (2001, 2004) and Stevens et al. (2004, 2005) from SCUBA 
observations of X-ray obscured AGN (spectroscopically classified 
as broad-line quasars) at 1 < z < 3, which appear to have higher 
sub-mm output than a matched sample of unabsorbed quasars. The 
emerging observational picture is that a significant fraction of the 
long-wavelength radiation has to originate in a strong starburst, 
thus suggesting that these objects are simultaneously building up 
their central SMBH, being obscured by large amount of dust and 
gas, and stars, the estimated duration of the obscured growth phase 
being ~ 15 per cent of the unobscured phase (Page et al. 2004; 
Stevens et al. 2005). If this picture were true, the unification model 
for AGN (Antonucci 1993), which explains the Type 1 (i.e., broad- 
line) vs. Type 2 (i.e., narrow-line) AGN dichotomy as due to geo- 
metric effects, would necessarily need to be considered in a broader 
context, where AGN evolution and feedback mechanisms are im- 
portant elements (e.g., Menci et al. 2008). 

To investigate the topics discussed above, we chose a sam- 
ple of 16 luminous obscured quasars at z ~l-2 detected in the 
2-10 keV band in the HELLAS 2XMM survey (see fl2~fl >, for which 
Spitzer follow-up observations have been performed (Pozzi et al. 
2007; Pozzi et al., in preparation). Approximately half of this sam- 
ple has a spectroscopic redshift and classification (Fiore et al. 2003; 
Maiolino et al. 2006; Cocchia et al. 2007); while IRAC and MIPS 
(24 fim) data are available for the entire sample, Spitzer observa- 
tions at 70 iim and 160 vim have been carried out only for eight 
sources with optical spectra. Finally, sub-mm observations have 
been performed for four of these sources only. 

Among the sample of 16 X-ray selected sources, 
H2XMMJ 003357.2-120038 (RA=00 33 57.2; DEC=-12 00 38.3; 
hereafter GD 158.19), a Type 2, radio quiet (Si. 4 GHz < 0.13 mjy 
at the 3o level) quasar at z= 1.957, because of its broad-band 
coverage (up to 850 (j.m) with good-quality photometric data, 
represents an ideal target to investigate the presence of both 
obscured AGN growth and star-forming activity at high redshift 
through a multi-band approach, and to estimate the bolometric 
energy output related to the accretion activity. In the following 
we will report on the multi-wavelength observational campaign 
conducted by our group over the last four years (§2) and on the 
modeling of the spectral energy distribution (SED) of this source 
(§3). The results will be discussed in §4 and summarized in §5. 

We note that GD 158_19 represents one of the few X-ray se- 
lected obscured AGN at high-redshift with such multi-wavelength 
characterization [in particular, with sub-mm detection; see also 
Mainieri et al. 2005 for a z=3.660 AGN in the Chandra Deep Field- 
South (CDF-S)]. 
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Figure 1. Optical spectrum of GD 158.19 smoothed with a boxcar of 16A; 
the main emission features are labeled (ESO, EFOSC2, 30-minute expo- 
sure). 



Hereafter we adopt Hq=70 km s Mpc~' in a A-cosmology 
with £1m=0.3 and O. A =0.7 (Spergel et al. 2003). Magnitudes are 
reported in the Vega system. 



2 MULTI- WAVELENGTH DATA REDUCTION AND 
ANALYSIS 

2. 1 GD 158.19 in the context of the HELLAS2XMM survey 

GD 158.19 was detected by the HELLAS2XMM survejQ, which 
comprises serendipitously detected sources in archival XMM- 
Newton observations over (originally) 3 deg 2 (Baldi et al. 2002), 
with 2-10 keV fluxes in the « 6x 10~ 15 -4x 10~ 13 erg cm~ 2 s~ ! 
range. The final source catalog of 232 X-ray sources over 1.4 deg 2 
(214 with optical counterpart down to R-band magnitude of ~ 25) 
is ~ 70 per cent spectroscopically identified (Cocchia et al. 2007) 
and comprises both Type 1 and Type 2 AGN, emission-line galax- 
ies, and early-type galaxies with indications of AGN activity from 
X-rays (XBONGs; see Comastri et al. 2002; Civano et al. 2007). In 
particular, albeit of low quality, the X-ray spectrum of GD 158.19 
was found consistent with being obscured by a rest-frame col- 
umn density of 7.3±,! j 1 5 7 x 10 22 cm~ 2 (Perola et al. 2004), once 
the photon index F was fixed to a "canonical" AGN value of 
1.9 (e.g., Piconcelli et al. 2005). The observed 2-10 keV flux of 
2.3xl0~ 14 erg cm~ 2 s _1 , corresponding to an intrinsic rest-frame 
2-10 keV luminosity (~ 6.3 x 10 44 erg s _1 ), coupled with the pres- 
ence of X-ray absorption, classifies GD 158.19 as an X-ray lumi- 
nous obscured quasar. Finally, we note that the X-ray luminosity 
is slightly higher (~ 20 per cent) than L* (at z = 2) from the La 
Franca et al. (2005) X-ray luminosity function. 

2.2 Optical spectroscopic classification 

Optical (R-band) photometry and long-slit spectroscopy (grism 
#13, covering a range of ss 3700-9300 A) of GD 158.19 were ob- 
tained in 2001 at the ESO 3.6m (EFOSC2 instrument), for a total 
exposure time of 10 minutes and 30 minutes, respectively; for de- 
tails, see Fiore et al. (2003). The object is relatively bright in the 
optical (R=21.8); all of the clearly identified permitted emission 

1 See the web page http://www.bo.astro.it/~hellas/sample.html 




lines (C IV and He II) have a full width at half maximum (FWHM) 
<1500 km s , In the noisy, red part of the spectrum (see Fig. [T) 
there are some indications of a broad Mg II emission line, which is 
not centered at the expected rest-frame wavelength of 2798A. The 
emission-line shift between the broad Mg II and the narrower C IV 
would be unusually large (~ 9000 km s _1 of redshift of the Mg II 
with respect to the emission features in the blue part of the spec- 
trum). This shift is three times as large as the maximum shift found 
by Richards et al. (2002) in a large sample of SDSS quasars. The 
current data suggest that either some spurious contribution from 
bad pixels redward of the Mg II line or residuals from background 
subtraction (or both) mimic a spurious broad emission line. Fur- 
thermore, also the redness of the continuum, the equivalent widths 
of the lines and their ratios, which are different from the values ob- 
tained from the composite Type 1 quasar spectrum obtained from 
the SDSS (Vanden Berk et al. 2001), are all suggestive of a Type 2 
AGN. We also note that GD 158_19 has different properties with re- 
spect to those of the X-ray obscured quasars at high redshifts found 
by Page et al. (2001, 2004; see also Stevens et al. 2005), which 
are characterized by broad-line optical spectra. Therefore, both the 
X-ray and optical properties place GD 158.19 in the long sought- 
after class of "genuine" Type 2 quasars (see also Comastri, Vignali 
& Brusa 2002). 



2.3 Spitzer data 

GD 158_19 was observed by Spitzer in 2006 with both IRAC and 
MIPS instruments in photometry mode, for a total integration time 
of 480 s (IRAC), 1400 s (MIPS at 24 urn), 3000 s (MIPS at 70nm), 
and 600 s (MIPS 160 |im). For the IRAC bands, we used the fi- 
nal post-calibrated data (PBCD) produced by the Spitzer Science 
Center (SSC) pipeline. For the typically more "problematic" MIPS 
bands, we started the analysis from the basic calibrated data (BCD) 
and we constructed our own mosaics using the SSC MOPEX soft- 
ware (Makovoz & Marleau 2005). At 24 iim, the input BCD for 
MOPEX were previously corrected using "ad-hoc" procedures to 



Table 1. GD 158.19: optical-sub-mm photometry. 



Observed 


Instr. 


Flux Density 


Band 




C"Jy) 


R 


ESO/3.6m 


5.6±0.8 


3.6 iim 


Spitzer/IRAC 


226±23 


4.5 \xm 


SpitzerflRAC 


387±39 


5.8 |om 


Spitzer/IRAC 


758±76 


8 |im 


Spitzer/IRAC 


1547±155 


24 nm 


Spitzer/MlPS 


5326±534 


70 \xm 


Spitzer/WBPS 


270()0±4500 


160 \xm 


Spitzer/MlPS 


37700±7800 


450 \im 


SCUBA 


< 94300 


850 (jm 


SCUBA 


8610±2060 



The upper limit to SCUBA 450 nm flux density is 
reported at the 3o" confidence level. 



remove the residual flat fielding, dependent on the scan mirror po- 
sition (see also Fadda et al. 2006; Pozzi et al. 2007). At 70 and 
160 p.m, we used the median high-pass filtered BCD (FBCD) to 
produce the mosaic map, as suggested for faint point-like sources 
(see the MIPS Data Handbook, Version 3.3.0). 

The flux densities in the IRAC and MIPS bands were mea- 
sured on the signal maps using aperture photometry. The chosen 
aperture radius for the IRAC bands was 2.45" and the adopted fac- 
tors for the aperture corrections are 1.21, 1.23, 1.38 and 1.58 (fol- 
lowing the IRAC Data Handbook, Version 3.0) at 3.6 iim, 4.5 \im, 
5.8 p.m, and 8 [im, respectively. 

The chosen aperture radius for the MIPS bands was 7" at 
24 iim and 16" at 70 iim and 160 nir0; me aperture flux densi- 
ties were corrected adopting a factor of 1.61, 2.07 and 4.1 in the 
three MIPS bands, respectively (see the MIPS Data Handbook). At 
70 p.m and 160 |im, these factors are well suited for sources with 
dust temperature of T~50-60 K. At these wavelengths, we were 

2 The PSF FWHM in the MIPS bands is ra 6" (24 yun), 18" (70 yim), and 
40" (160 iim); see MIPS Data Handbook. 



forced to use a small extraction radius (in particular at 160 [im, 
where the radius is comparable to half of the PSF FWHM), since 
at a distance of fa 50" from our target there is a source (undetected 
in the X-ray band) of comparable brightness in the IRAC bands 
but fa 35 and 20 per cent brighter at 70 [im and 160 ixm, respec- 
tively (Fig. [2}- The companion source would contaminate the target 
emission if larger apertures were adopted. Using an aperture diam- 
eter smaller than the PSF FWHM (as in the case of 160 iim data) 
means to rely on the PSF core, where flux non-linearities could 
play a significant effect due to the limited sampling of the PSF (see 
Gordon et al. 2007). On the other hand, we have tested how a PSF 
fitting would not provide reliable results, probably due to the small 
field-of-view of the mosaics (fa 2-3 arcmin; see also Stansberry et 
al. 2007, where the latest 160 |J.m flux density absolute calibration 
is discussed). 

Table [Tj reports the target flux densities provided by Spitzer. 
To compute the uncertainties, we added in quadrature the noise map 
and the systematic uncertainty, assumed conservatively to be 10 per 
cent in the IRAC and MIPS 24 iim bands, and 15 per cent at 70 iim 
and 160 \im (see IRAC and MIPS Data Handbooks). While the 
systematic uncertainties dominate in the IRAC and MIPS 24 p_m 
bands, at longer wavelengths also the noise map is an important 
factor of uncertainty, approaching the level of the systematic un- 
certainty at 160 |i.m. 

We note that, given our flux density measurements in the 
IRAC bands, GD 158_19 would have been recognized as an AGN 
in both the Lacy et al. (2004) and Stern et al. (2005) diagnostic 
diagrams. Also the companion source would lie in the AGN "lo- 
cus" in these diagrams. Given the source properties in the mid-IR 
and far-IR, it appears plausible that GD 158_19 and its companion 
are both characterized by AGN and starburst emission. The X-ray 
non-detection of the companion source may indicate either weak 
X-ray emission or obscuration, although the limited exposure of 
the XMM-Newton observation (fa 21.6 ks and 13.8 ks in the MOS 1 
and MOS2, respectively) prevents from further investigation. 

Using the 24 [im galaxy counts in the SWIRE fields from 
Shupe et al. (2008), we computed the integral number density of 
24 \im sources with flux density above fa 5.3 mJy (the flux densi- 
ties are 5.33 mJy and 5.20 mJy for GD 158_19 and the companion, 
respectively), which is ~ 25 deg~ 2 (a similar number density is 
found in Papovich et al. 2004); this number translates into 0.015 
expected sources in a circular region with radius of 50" (the dis- 
tance between the two sources). The Poisson probability of obtain- 
ing one or more sources when 0.015 are expected is relatively low 
(?a 1.49 x 10~ 2 ); if placed at the same redshift, they could repre- 
sent the density peaks in an over-dense region, as found in the fields 
around some high-redshift radio galaxies with associated powerful 
sub-mm emission (e.g., Stevens et al. 2003; Smail et al. 2003). 

2.4 SCUBA data 

The SCUBA observations of GD 158.19 were split in two nights 
over October-November 2004 and carried out in photometry mode, 
in which the source is placed on the central bolometer of the array, 
with the wide 850:450 filter set. The sky opacity was monitored 
on a quasi-continuous basis with the JCMT water vapour radiome- 
ter and cross-checked with skydips several times per night. Uranus 
was chosen as flux-density calibrator; the estimated calibration un- 
certainties on the most sensitive 850 \im data are fa 10%. Data 
were reduced with the Starlink SURF package in a standard manner 
(e.g., Holland et al. 1999). We reduced data from the two nights 
separately and then calculated a weighted mean flux density and 



error. The source was detected at the 4.2a level at 850 iim, while 
at 450 iim only a loose upper limit was obtained (reported at the 
3(7 level in Table[T). Although the Spitzer 160 [im data reveal the 
presence of a bright infrared "companion" source, its distance from 
GD 158_19 is such that it should not give any significant flux contri- 
bution to our target at 850 |J.m given the SCUBA beam size (fa 15") 
at this wavelength. 



3 DATA MODELING 

The optical-to-sub-mm data of GD 158.19 have been modeled us- 
ing the code developed by Fritz et al. (2006, hereafter F06). This 
code assumes a continuous dust distribution (around the central 
source) consisting of silicate and graphite grains, confined in a 
toroidal shape (Pier & Rrolik 1992), as opposed to the "clumpy" 
tori models mostly developed over recent years (e.g., Nenkova, 
Ivezic & Elitzur 2002; Honig, Beckert & Ohnaka 2006; Nenkova 
et al. 2008a,b; Schartmann et al. 2008). Notwithstanding this prob- 
ably simplified underlying modeling, the fundamental goal of the 
present work consists in estimating the physical parameters related 
to the AGN (torus) emission, peaking in the mid-IR domain. We 
adopted the multi-component fit approach (e.g., Hatziminaoglou et 
al. 2008); in the following, we resume its main characteristics. 

We use three distinct components that, in general, dominate 
different parts of a galaxy SED: 

(i) stellar emission (UV-near-IR); 

(ii) emission due to hot dust heated by AGN emission (mid-IR); 

(iii) emission due to warm/cold dust heated by star formation 
(starburst component, mid- and far-IR). 

(i) To model the stellar component, we use a set of Simple 
Stellar Population (SSP) spectra of solar metallicity and ages rang- 
ing from fa 1 Myr to 1.7 Gyr; the latter is the time elapsed between 
z=4 (the redshift assumed for the stars to form) and z = 1.957 in the 
adopted cosmology. A Salpeter (1955) initial mass function (IMF) 
has been assumed. These are weighted by a Schmidt-like law of 
star formation: 

^ W = ^exp(!_l) (1 ) 

where Tq is the age of the oldest stellar population and t sets the 
duration of the initial burst (normalized to Tq). A common value 
of extinction is applied to stars of all ages, and the extinction law 
of the Galaxy has been adopted (Ky =3.1; Cardelli, Clayton & 
Mathis 1989). 

(ii) Type 2 AGN emission, dominating in the mid-IR domain 
(from a rest-frame wavelength of fa 2 up to fa 30-50 (J.m) and 
produced in a dusty region surrounding the SMBH, is modeled us- 
ing an extended version of the model grid described by F06. 

(iii) The emission from warm/cold dust, peaking in the far- 
IR (at fa 100 |i.m rest frame) and being characterized by features 
in the mid-IR due to PAH molecules, has been accounted for by 
using semi-empirical models of known and well studied starbursts 
(Arp 220, M 82, M 83, NGC 1482, NGC 4102, NGC 5253 and 
NGC 7714). 

The fit to the observed photometry is carried out at first on 
the UV-optical-near-IR region (over the fa 0.3 — 5 ixm rest-frame 
interval), by means of a pure stellar component. Then all the star- 
burst templates are normalized to match the far-IR data (longward 
of fa 25 \im rest frame). At this point, the best AGN model is 
searched, among all of those available from the aforementioned 
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Figure 3. Observed-frame SED of GD 158.19 from the R-band to the 
SCUBA 850 |am data points. The upper limit refers to the source non- 
detection at 450 [Mil (SCUBA data). The "global" SED (solid line), re- 
produced as the summed contribution of a stellar component (dotted line), 
a torus component (long-dashed line), and a starburst component (short- 
dashed line), represents the best-fitting according to F06 modeling. 

model grid, by varying the torus parameters (e.g., effective angle, 
viewing angle, dust distribution, etc.). Finally, the normalizations 
of the three components are iteratively changed until a good global 
solution is achieved; at this stage, also the possibility that one of 
these three components may not be necessary is taken into account. 
The goodness of the fit is evaluated by means of a "merit function" 
which may be intended as a modified reduced % 2 , where the num- 
ber of degrees of freedom is replaced by the number of photometric 
points used in the SED fitting, following the equation 

MF = Z[-^- L ) l N o (2) 

where No is the number of observed datapoints, M,- and O, are the 
model and observed flux densities, respectively, of the i — th photo- 
metric band, and <J, is the corresponding observed error. Since the 
observed SED at study includes only photometric data and, further- 
more, is not fully sampled, the problem of finding a best-fit model 
is intrinsically prone to a certain degree of degeneracy, i.e. to the 
non-uniqueness of a solution. 

As a posteriori check, we estimated what cx ox (defined as 
the slope of a hypothetical power-law connecting the rest-frame 

2500 A to 2 keV (oc ox = '° g | /2tcV ^ 25QoA) , ) is needed to reproduce 

log(v, kcV/V 2500 A ) * 

the source rest-frame, de-absorbed 2-10 keV luminosity. This will 
allow us to provide a comparison with recent literature results based 
on optically selected broad-line quasars (e.g., Vignali, Brandt & 
Schneider 2003; Steffen et al. 2006). 



4 RESULTS FROM THE SED FITTING AND 
DISCUSSION 

The output of the SED fitting procedure described above consists 
of a series of models with different MF, each having its own set of 
parameters describing the stellar, AGN and starburst components. 
In the following, we will describe the best-fitting model parameters 
related to the SED of GD 158_19. 

The best fitting has a MF of ~ 2.2; most of the discrepancy 



is due to the modeling over the wavelength range 24—70 \im (see 
Fig. [3}, where the model SED is a factor ~ 2 below the photometric 
points. However, we note that apparently unacceptable MFs repre- 
sent a common problem of most of the SED fittings (see, e.g., Fig. 2 
of Hatziminaoglou et al. 2008), due also to the combination of pho- 
tometric measurements with largely different accuracies (see, for a 
comprehensive description of this issue, Gruppioni et al. 2008). 

The AGN component in the best-fitting "global" SED model 
(solid line in Fig. |3) consists of reprocessed emission by an al- 
most edge-on torus (~ 10 degrees with respect to the equatorial 
plane) - as expected from the optical and X-ray classifications of 
GD 158_19 as a Type 2 quasar - with a full covering angl^l of 
~ 140 degrees (corresponding to a covering factor of ss 90 per cent) 
and an optical depth (at 9.7 (im) of t(9.7) = 1.0, corresponding to 
a column density of ~ 9xl0 22 cm~ 2 for a Galactic dust-to-gas 
ratio. This value is similar to the one derived (though with large 
uncertainties) from the X-ray spectrum. Given the source X-ray lu- 
minosity, the corresponding optical-to-X-ray spectral slope would 
be oc ox =- 1.55 ^-1.60. 

The AGN contribution to the rest-frame 1-1000 \tm emis- 
sion is 54 per cent; the remaining emission is mostly ascribed to 
a prominent starburst emission, best reproduced by an Arp 220 
template (Berta et al. 2003; short-dashed line in Fig. [3}, made by 
combining GRASIL (Silva et al. 1998) models and the 5-20 \im 
ISOCAM-CVF spectrum of the source (Charmandaris et al. 1999). 
Despite this template under-estimates the observed flux density at 
70 [J.m (close to the peak of the torus emission), we note that no 
other starburst template is able to reproduce the far-IR and sub-mm 
emission of GD 158_19 better than Arp 220. 

From the infrared luminosity - once the AGN contribution has 
been removed - and adopting the Kennicutt (1998) relation, we ob- 
tain a SFR of ~ 1500 Mq yr~' . For comparison and further check, 
the recent SFR was also computed using Eq. (1), i.e., integrating the 
law of star formation through fitting the rest-frame optical/near-IR 
datapoints with the stellar component (dotted line in Fig. [3}, char- 
acterized by E(B— K)~0.49. The derived SFR is consistent with 
that obtained from the Kennicutt (1998) relation. 

The SFR for GD 158.19 appears to be consistent with the val- 
ues derived by Stevens et al. (2005) for a sample of broad-line, 
X-ray obscured quasars detected with SCUBA over the redshift 
range z~l-3 and with the sub-mm detected quasars at z ~ 2 from 
Coppin et al. (2008). For comparison, we note that the SFR de- 
rived by Mainieri et al. (2005) for an X-ray selected Type 2 quasar 
at z=3.66 is a factor ~ 2.2-2.5 lower. At zero-th order, there are 
significant indications that GD 158.19 is experiencing a phase of 
strong star formation and AGN activity, similarly to other quasars 
at z » 2. 

No other SED fitting solution using the three-component mod- 
eling provides similar results in terms of MF and overall represen- 
tation of the X-ray properties of GD 158.19. While in all of the 
SED fittings the full covering angle of the torus is 140 degrees, so- 
lutions with x(9.7) < 1.0 (corresponding to higher inclination an- 
gles for the line of sight with respect to the equatorial plane of 
the torus) would significantly under-estimate the observed column 
density. Similarly, a solution requiring a flatter oc ox (—1 .30, close to 
the average oc ox found in X-ray selected broad-line quasars in the 
XMM-COSMOS survey; Lusso et al., in preparation) would imply 
a worse MF (ss 4.4), mostly due to its inability in reproducing the 
source emission over the wavelength range 5-8 vim (where the stel- 

3 We note that in F06 this angle was referred to as full opening angle. 



lar and AGN components are approximately equal); in this case, the 
parameters derived for the AGN and starburst emission are similar 
to those presented above. We also note that solutions with only the 
AGN plus starburst emission can be discarded on the basis of the 
low optical depth (t(9.7)=0.09, i.e., a factor of ~ 10 lower column 
density than actually measured); in this case, a steeper oc ox , close 
to ~ —1.8, would be required to reproduce the observed SED, and 
a factor ~ 4 higher bolometric luminosity than that obtained with 
the best-fitting solution (see below) would be implied. 

Despite the fact that t(9.7) along the line of sight is equal 
to 1, the 9.7 (j.m silicate feature is seen in emission in the best- 
fitting model, while one would expect it in absorption. This is due 
to the combination of two effects: first, having a lower sublimation 
temperature with respect to graphite, the silicate grains can only 
survive at larger distances from the inner radius of the torus. Thus 
the optical depth, calculated from this "hot-silicate" zone towards 
the observer, can be, in some cases, lower than unity (because not 
integrated throughout the entire thickness of the torus), thus mak- 
ing possible the observation of this feature in emission. Second, in 
models where the density is decreasing at increasing height over the 
equatorial plane, the emission feature comes from the lower density 
zones, where the hot silicates are living in an optically thinner en- 
vironment. 

According to recent results, a limited fraction of optically lu- 
minous quasars at z ~ 2 are detected in the sub-mm and millimeter 
(mm) wavelengths (~ 25 per cent; Omont et al. 2003), suggesting 
that the quasar and sub-mm phases do not overlap significantly, 
probably because of the different lifetimes of the two phases. 
Within this scenario, a quasar detected in the sub-mm/mm is a good 
candidate for being caught in the so-called "transition phase" from 
a growing, sub-mm galaxy to an unobscured quasar. In such a phase 
(whose duration is estimated of the order of ~ 15 per cent of the un- 
obscured phase; Stevens et al. 2005), powerful obscured emission 
from the active nucleus, coupled with a prominent starburst com- 
ponent, is expected and actually observed (e.g., Page et al. 2001, 
2004; Stevens et al. 2005; see also Aravena et al. 2008). 

We estimate an AGN bolometric luminosity, integrated over 
the optical, UV and X-ray band (extrapolated up to 500 keV assum- 
ing a power law with photon index T=1.8 and a cut-off at 300 keV), 
of ?s4.3xl0 46 ergs" 1 . 

Since the black hole mass of GD 158.19 cannot be derived 
using "standard" methods (e.g., reverberation, scaling relationship, 
etc.), we decided to estimate it using the MBH-Mt, u i ge relation de- 
rived locally by Marconi & Hunt (2003) from a sample of ~ 30 
galaxies with a "secure" black hole mass measurement (hence, no 
evolution of the relation over cosmic time is assumed). The stel- 
lar mass, obtained directly from the best-fitting SED modeling, 
is ~ 1.5 x 10 12 M Q , close to the upper boundary for the stellar 
masses reported by Dye et al. (2008; see their Fig. 9) for a sample 
of SCUBA galaxies at z ~ 2. The assumption of a Chabrier (2003) 
IMF, which has been suggested to be more representative of the 
IMF at high redshifts (e.g., Renzini 2006), would imply a factor 
~ 1.7 lower stellar masses, without significant changes in the over- 
all SED fitting. The estimated black hole mass is fa 1.9 x 10 9 M G , 
implying an Eddington ratio of ~ 0.19, higher than the values 
recently found for a sample of HELLAS2XMM selected Type 2 
quasars at high redshift (~ 0.05 on average; Pozzi et al. 2007; Pozzi 
et al., in preparation; for further details on the source selection, see 
also £D- 

If we assume that the star-formation episode lasts, at the rate 
estimated in this paper, for a time interval t, while the black hole 
grows at the current rate for a period / qso = T| x t (where r| de- 



pends on the presence of available fueling for accretion, hence on 
the history of mergers and dynamical instabilities of GD 158.19), 
we derive that the system would move ~ 0.3 dex above (below) the 
Marconi & Hunt (2003) relation over a timescale of ~ 6 x 10 s yr 
(~ 4 x 10 9 yr) if T|=1.0 (0.1) is assumed. In other words, the per- 
sistence of the system along the relation depends on the relative 
duration of the quasar and starburst; in the case of T|=1.0, that rep- 
resents the situation in which their phases last for a comparable 
amount of time (at the rates estimated in our work), it would take 
less than 1 Gyr to move the system away from the relation. 

An alternative, and admittedly speculative, possibility is that 
the cold gas and dust reservoir, which sustain the star formation 
at the estimated rate of ~ 1500 M@ yr -1 , is also responsible for 
the obscuration of the nuclear source. The observed X-ray absorp- 
tion can be easily reproduced for a relatively large range of the 
Lamastra, Perola & Matt (2006) molecular disk parameters, pro- 
vided that the surface density exceeds ~ 8 M© pc~ 2 . We note that 
the molecular gas has to be close to the nuclear region, otherwise 
the narrow-line region would be hardly visible. 



5 SUMMARY OF THE RESULTS 

We have presented the optical, infrared and sub-mm observations 
of a luminous, X-ray obscured Type 2 quasar at z= 1.957 detected in 
the 2-10 keV band within the HELLAS2XMM survey. By fitting 
the photometric points with a three-component fit approach (stellar, 
AGN and starburst emission; see F06 for details), we were able to 
estimate some fundamental parameters related to the physical and 
geometrical properties of the torus and starburst components. 

• The AGN component, related to emission reprocessed by dust 
in the torus, is dominant in the mid-IR, peaking at ~ 20-30 p.m 
(rest frame). The obscuring matter is likely distributed in an al- 
most edge-on structure (~ 10 degrees with respect to the equatorial 
plane), with a full covering angle of ~ 140 degrees. The derived 
optical depth [x(9.7) = 1.0], corresponding to a column density of 
~ 9x 10 22 cm~ 2 , is in agreement with the absorption derived by 
the X-ray spectral fitting. 

• The AGN contribution to the rest-frame 1-1000 p.m emission is 
of ~ 54 per cent, the remaining being ascribed to a starburst emis- 
sion (best reproduced by an Arp 220 template) peaking at longer 
wavelengths with respect to the AGN emission. The derived SFR, 

1500 M yr~' , places GD 158.19 as a breeding ground for con- 
siderable star formation at z ~ 2. 

• The picture emerging from the results of the SED fitting is that 
of an object caught in a phase where intense, dust-enshrouded star 
formation is coeval to obscured AGN growth, which however does 
not prevent the active nucleus from revealing its presence in the 
optical and X-ray bands. If current evolutionary scenarios (e.g., 
Hopkins et al. 2006; Menci et al. 2008) were correct, we could 
expect that the black hole in GD 158.19 is in the process of be- 
coming a "classic", luminous and unobscured quasar, capable of 
quenching the star formation. There are indications that such phase 
is shorter for the most luminous AGN (e.g., Adelberger & Steidel 
2005; Menci et al. 2008). GD 158.19 seems to have properties in 
common with the sub-mm detected and X-ray obscured quasars at 
z~l-3 found recently in other surveys (Page et al. 2004; Stevens et 
al. 2005). 
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